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RADIOMETRIC MEASUREMENTS OF THE EARTH'S INFRARED HORIZON 
FROM THE X-I5 IN THREE SPECTRAL INTERVALS 

By Antony Jalink, Jr., Richard E. Davis, 
and Dwayne E. Hinton 
Langley Research Center 

SUMMARY 

Experimental data of the earth's horizon in three spectral intervals, 0.8 to 2.8 pm, 

10 to 14 pm, and 14 to 20 pm, have been obtained by using the X-15 research airplane. 

The resolution (field of view) of the data at the tangent point is approximately 2 kilome- 
ters, and the spatial position of the horizon profile is known to within 4 kilometers. 

The 14 to 20 pm interval was shown to be the best spectral interval for horizon 
sensing because it exhibits the least variance at adequate radiant intensity. The 0.8 to 
2.8 pm and 10 to 14 pm spectral intervals do not appear suitable for accurate horizon 
sensing because they are sensitive to low-altitude meteorological changes. 

The experimental data obtained verified, within experimental accuracy, the theory 
used to predict the far infrared radiant energy due to carbon dioxide. The use of theory 
for spectral regions where high- altitude water-vapor data are required is questionable 
because no mixing- ratio information is available. The single- scattering theory cannot 
predict accurately the radiant energy in the near infrared region {<3.0 pm). 

INTRODUCTION 

Conventional flight- vehicle attitude instruments use gravity sensors to determine 
the local vertical. In orbital vehicles, which are in essentially a "zero g" condition, 
such instillments become useless. As an alternative, when the space vehicle is in the 
vicinity of the earth, the apparent horizon formed either by thermal emission or by 
scattered and reflected sunlight from the earth and its atmosphere can be sensed against 
the effective temperature of interstellar space (3.1° K) to determine vehicle attitude. 
Ideally, the radiation profile sensed by a horizon scanner when it scans from space onto 
the earth should be well defined; that is, it should be steep and invariant. However, the 
radiation gradient is affected by meteorological conditions, solar ir radiance, location on 
earth, and the spectral interval employed by the sensor. Knowledge of the effect of these 
factors on the infrared characteristics of the earth's radiance profile determines, to a 
great extent, the positional accuracy and reliability attainable by horizon- sensing methods. 


In order to obtain experimental measurements of the radiation gradient, a program using 
the X-15 research airplane as a test vehicle (ref. 1) was instituted by the Langley 
Research Center. A radiometer carried in the tail-cone box of the X-15 was used to 
view the horizon during the peak- altitude portion of the X-15 flight. 

This report presents the horizon- radiance measurements obtained from three 
flights. The spectral band and data of each flight is as follows: 


Flight number 

Spectral band 

Date 

Flight 1 

0.8 to 2.8 jam 

July 8, 1964 

Flight 2 

10 to 14 pm 

May 28, 1965 

Flight 3 

14 to 20 pm 

June 29, 1965 


Meteorological measurements were obtained for each flight to allow construction of 
calculated radiance profiles based on existing radiance theory for comparison with the 
measured profiles. This report contains a description of the flight instrumentation 
employed for the X-15 missions, which includes the design and calibration of the radiom- 
eter, the description of the X-15 instrumentation used to obtain the measurements, and 
the method used to calculate the radiance profile relative to the earth's surface. A dis- 
cussion of the instrumentation errors and the expected accuracy of the measurement sys- 
tem is given. 
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SYMBOLS 

area, centimeters^ (cm^) 
equivalent noise bandwidth, hertz (Hz) 
diameter of objective, centimeters (em) 
detectivity, centimeter-hertz 1 ^ 2 /watt (cm-Hz^/w) 

water-vapor pressure, millibars (rabar) (l millibar = 1 hectonewton/meter2) 

focal length of collimator, centimeters (cm) 

focal length of radiometer optics, centimeters (cm) 

frequency, hertz (Hz) 


2 



corner frequency, hertz (Hz) 


fc 

G(jaj) transfer function 

g acceleration due to gravity, meters/ second 2 (m/ sec 2 ) 

H irradiance, watts/meter^ (W/m2) 

h tangent height (see fig. 1 ), kilometers (km) 

h x _i 5 altitude of X-15, kilometers (km) 

K midband gain 

Kj,Kg constants 

l distance, centimeters (cm) 

Mg, M 3 , M 4 matrices 

N number of samples 

n number of samples 

P\ spectral radiant power, watts/micrometer (W/ /j.m) 

p total atmospheric pressure, millibars (mbar) 

R radiant emittance, watts/meter2 (w/m2) 

Rx spectral radiant emittance, watts/meter^-micrometer (W/m^-fim) 

R^ spectral radiance, watts/meter^-steradian-micrometer (W/m^-sr-pm) 

r radius of earth, kilometers (km) 

T temperature, degrees Kelvin (°K) 

t time, seconds (sec) 
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times of evaluation, seconds (sec) 
initial voltage, volts (V) 
noise voltage, volts (V) 
signal voltage, volts (V) 

water-vapor mixing ratio, grams/kilogram (g/kg) 

difference between average value and value of sample N 

local earth axes 

aircraft axes 

radiometer axes 

space-fixed axes 

scan mirror angle, degrees (deg) 

vertical mounting angle, degrees (deg) 

nadir angle, degrees (deg) 

optical efficiency of system 

pitch angle, degrees (deg) 

error in pitch angle due to platform calibration and recording paper shrink- 
age, degrees (deg) 

minimum angle of resolution, degrees (deg) 
wavelength, micrometers (/im) 
standard deviation 
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a h standard deviation in tangent height, kilometers (km) 

^X-IS standard deviation in X-15 altitude, kilometers (km) 

a a standard deviation in scan mirror angle, degrees (deg) 

<7^ standard deviation in vertical mounting angle, degrees (deg) 

Og standard deviation in pitch angle, degrees (deg) 

0 g c standard deviation in error in pitch angle, degrees (deg) 

cr^ standard deviation in roll angle, degrees (deg) 

t time constant, seconds (sec) 

<#> roll angle, degrees (deg) 

\p yaw angle, degrees (deg) 

horizontal mounting angle, degrees (deg) 
io solid angle, steradians (sr) 

a> v field of view, steradians (sr) 

EXPERIMENTAL SYSTEM DESCRIPTION 

The horizon profile can be represented by a plot of radiance as a function of tangent 
height, where tangent height is the minimum altitude of the radiometer line of sight. Tan- 
gent height can be computed from the position of the observer and the nadir angle to the 
line of sight of the radiometer, as shown in figure 1 and derived in appendix A. The nadir 
angle, in turn, is found from the instantaneous attitude of the aircraft and the angle of the 
mirror scan. If the calculations are performed in all degrees of freedom, the latitude 
and longitude of the tangent point are also obtained. 

A systems diagram of the X-15 horizon definition experiment is shown in figure 2. 
The dashed-line boxes indicate available data sources routine to X-15 flights. The solid- 
line boxes show equipment and calculations provided especially for this experiment. The 
radiometric output is combined with the calibration data obtained for the instrument so 
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that the observed radiance is known. From the X-15 stable platform, aircraft pitch, roll, 
and yaw data are obtained. These data are combined with the radiometer scan mirror 
angle and aircraft position from ground tracking to obtain the direction and position of the 
instrument line of sight so that the observed radiance data can be fixed relative to a spe- 
cific point on earth. A set of data points is extracted from the flight records every 
5 msec, and the calculations are performed with the aid of an IBM 7094 computer. 
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Figure 2*- A systems diagram of horizon definition experiment. 

RADIOMETER DESIGN 

The radiometer is mounted in the tail-cone box, immediately behind the vertical 
stabilizer fin of the X-15 rocket airplane, and looks back out of the rear of the aircraft, 
as shown in figure 3. Because the tail-cone box is external to the airplane and is not 
heated, the radiometer is subject to a coldsoak near -50° C prior to launch when the X-15 
is carried to an altitude of 14 km by a B-52 aircraft. All structural parts of the radiom- 
eter are invar to minimize defocusing of the optics due to the extreme temperature range. 
Also, the X-15 rocket engine operation results in a lOg vibration environment which the 
instrument must survive. However, during the data period which occurs at peak altitude 
the engine is shut down and interference from engine vibration and exhaust is thus elimi- 
nated. A 16-mm motion-picture camera mounted beside the radiometer observed the 
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Figure 3.- Radiometer mounted on X-15. L -68-861 


presence of clouds and haze during the data-taking period. Postflight examination and 
calibration of the radiometer and the film from the camera indicated that combustion 
products from the rocket exhaust did not measurably contaminate the optics. 

A schematic drawing and a photograph of the radiometer are shown in figures 4 
and 5, respectively. The radiometer uses a motor-driven scan mirror which sweeps the 
field of view through a 30° optical scan at a speed of 20° per sec. A precision potentiom- 
eter is used to record the mirror position. The scan mirror reflects the energy onto the 
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Figure 4.- Schematic drawing of X-I5 radiometer. Field of view, 0.13° by 0,B°; 
diameter of objective, 12.7 cm. 


objective, a 12. 7- cm- diameter paraboloidal mirror with an f-number of 1, which focuses 
the energy onto the detector. The field of view which is determined by the size of the 
detector, was measured to be 0.13° by 0.13° at the half- response points. This detector 
size results in a radiometric resolution of 2 km at the horizon tangent point. Because 
the detectors are affected by changes in the ambient temperature, they are thermostati- 
cally controlled within a preset range. The radiation also passes through an optical band- 
pass filter used to select the spectral region to be investigated. The filter is not heated, 
because the shift in wavelength due to the expected ambient- temperature range is, at 
most, 0.8 percent which is thought to be negligible. The optical filter, detector, and elec- 
tronics are somewhat different for each spectral interval. The radiometer characteris- 
tics for each flight are listed in table I and the spectral-response curves are shown in 
figure 6. A detailed discussion of the characteristics of each detection system and the 
methods used to calibrate them are given in appendix B. 
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Figure 5,- Photograph o# X“I3 radiometer- L-68-862 


METEOROLOGICAL INFORMATION FOR THE FLIGHT EXPERIMENT 


In order to evaluate the horizon- radiance theory and its ability to predict radiance 
profiles, it is necessary to obtain information concerning the meteorological conditions 
existing at the time of each flight. With this information, theoretical radiance profiles 
may be calculated and compared with the measured radiance profiles. 

Many authors (refs. 2, 3, and 4) have formulated procedures to obtain atmospheric 
radiance profiles in the far infrared region. In the 10 to 14 pm window interval, most of 
the radiance is received from opaque bodies, such as the earth, and clouds, and is only 
slightly attenuated by the atmosphere. In the 14 to 20 pm interval, the carbon dioxide 
and water vapor (%€>) within the atmosphere provide thermal radiation which is 
subsequently attenuated by the atmosphere between the radiating element and the detector. 
In order to calculate the radiance received, it is necessary to know the amount of absorber 
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Wavelength* \ita 


Wavelength, urn 


Wavelength p 14m 


(a) O.B to 2,8 tim. 


lb) 10 to 14 pm. 


<c) 14 to 20 pm. 


Figure 6,' Radiometer spectral- response curves for each flight 


and the temperature and pressure distributions along the line of sight. Because a spheri- 
cally stratified atmosphere is assumed in the calculations for mathematical simplicity, 
the construction of a model atmosphere consists of the specification of temperature and 
mixing ratio as functions of pressure altitude. (See ref. 3.) The model atmospheres 
derived for this purpose are listed in tables IX, III, and IV, The methods used to derive 
these atmospheres and the meteorological measurements employed are described in 
appendix C. 

The theoretical infrared horizon profiles for the 10 to 14 jam and 14 to 20 jam 
spectral intervals (corresponding to flights 2 and 3, respectively) were computed by 
J. Alishouse of the National Environmental Satellite Center of the Environmental Science 
Services Administration (ESSA) by using the ESSA infrared radiance calculation program 
(ref. 3). The profile for the 0.8 to 2.8 jam spectral interval was obtained from 
reference 5. 
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FLIGHT-TEST RESULTS 


Several of the radiance profiles which were measured in the three spectral intervals 
are listed in tables V to VII. The radiance profiles, together with the respective cloud 
charts, are shown in figures 7 to 12. The measured radiance profiles are compared with 
theoretical profiles based on the meteorological conditions existing at the time of each 
flight. The radiance values of the horizon profiles presented in this paper are given in 
terms of effective radiance. To compare the actual values of radiance from the horizon 
with those shown, the actual values must be multiplied by the radiometer optical filter 
response. An estimate of the accuracy of the measurements can be obtained from 
appendix D. 


Flight 1 

Flight 1, which was conducted on July 8, 1964, obtained horizon- profile measure- 
ments in the 0,8 to 2.8 |im spectral interval. The primary radiant energy in this wave- 
length interval is reflected and scattered solar energy, as shown by curve A in figure 13. 
The scattering particles include air molecules and atmospheric aerosols, such as dust, 
smoke, salts, and water. The amount of scattered energy received by the sensor depends 
strongly on the angle between sensor line of sight and the direction of solar illumination, 
as well as on the wavelength of the radiation, and size, type, and density of the particles. 

Several of the radiance profiles measured during flight 1 are given in table V and 
figure 7, The solid- line curve in figure 7 is a smoothed profile obtained by averaging the 
measured radiance of all horizon crossings at each tangent height. The cross-hatched 
area about this curve represents the envelope containing all data, and the dotted curve is 
typical of the measured profiles obtained. 

Figure 7 shows that sunlight scattered by the atmosphere causes a smooth increase 
in the energy received by the radiometer as its line of sight sweeps from a tangent height 
of 50 to 20 km. A steep increase in radiance occurs slightly below a tangent height of 
20 km and leads to a plateau which extends to a tangent height of approximately 12 km. 
This feature is apparent in the curve representing a typical horizon profile but is subdued 
in the average curve because of variation in altitudes between individual profiles. This 
energy can be explained by solar radiation scattered by one or several aerosol constitu- 
ents which surround the earth at altitudes of approximately 14 km and 19 km, and which 
are invisible from the earth's surface because they have little optical thickness in the 
vertical direction (ref, 6). 

Below a tangent height of 10 km the profile is dominated by energy from cirrus 
clouds which covered the area of earth where the profiles were obtained. Successive 
clouds can be noted at tangent altitudes of 10 km and -1 km for the typical horizon 
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crossing shown in figure 7. The clouds as detected by the radiometer are consistent with 
the conditions shown by the 16 -mm motion-picture camera mounted next to the radiom- 
eter and by the weather patterns existing that day as shown in figure 8. The dashed-line 
box in figure 8 shows the area where the horizon crossings were obtained. The area of 
view was partly clear but high cirrus clouds covered much of the sky to the east. The 
maximum measured effective radiance shown in figure 7 is 52 W/m^-sr. The solar 
energy which the radiometer would measure when looking at clouds with an albedo of 1 

is 94,5 W/m^-sr, so that the measured effective radiance corresponds to an albedo of 
52 

~~r = 0.55 which compares favorably with available data. 

94.5 

Figure 7 also shows the theoretical horizon profile for primary Rayleigh scattering 
{ref. 5). This profile applies to the northern horizon toward which the radiometer was 
pointing. The measured data indicate that more solar energy is scattered toward the 
observer than the theoretical curve predicts. Because calculations were limited to pri- 
mary molecular scattering, the higher energy reaching the radiometer can be accounted 
for by multiple scattering and the presence of aerosols. 

The results of tests in the near infrared region indicate that the 0.8 to 2.8 pm inter- 
val is not a good region for horizon sensors. The radiance above a tangent height of 20 km 
results from scattering. The gradient in this region is shallow and, as noted in the mea- 
sured results, it is variable because of multiple scattering, varying aerosol content, and 
small changes in sun angle. Between 15 km and 20 km, the gradient is steeper but the 
variance in the boundaries of the data increases. Below a tangent height of 15 km the 
gradient is very steep, but clouds are extremely variable; therefore, this region does not 
possess the qualities of a steep invariant radiance profile. 

Flight 2 

Figure 13 shows that beyond 4.6 pm a horizon sensor in space would receive radia- 
tion mainly from the earth and its atmosphere. Therefore, a horizon sensor which oper- 
ates in a spectral interval in the far infrared region can sense the horizon on the dark 
(night) side as well as the sunlit (day) side of the earth and is virtually independent of the 
sun's position. 

For a far infrared sensor which includes no, or at most weak, atmospheric absorp- 
tion bands within its spectral region, the sensor receives its radiant energy from opaque 
bodies, such as the earth itself, or from clouds that may screen it. The amount of radi- 
ance received by such a sensor from the earth for clear conditions {curve B of fig. 13) 
differs greatly from that received from clouds for overcast conditions (curve D of fig. 13) 
and, therefore, a strong dependence on meteorological conditions is implied. The spectral 
interval from 10 to 14 pm is an example of a window region where atmospheric absorption 
is low. 
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Figure fL- Mephanalysis (cloud chart} for GJ to L8 spectral interval Flight 1; C0:00 GMT, July 9 , 1964. 






Flight 2 was made on May 28, 1965, and measured the earth's horizon profile in this 
spectral interval. The data are shown in figure 9 where the solid curve is the averaged 
radiance profile as previously described. Other measured radiance profiles, which are 
all within the cross-hatched area, are given in table VI. The dashed curve is a theoreti- 
cal profile calculated by ESSA for the 10 to 14 pm radiometer spectral response shown 
in figure 6(b) and from the meteorological data for flight 2 . The measurements were 
obtained in an almost clear and cloudless area, as shown in figure 10, and represent the 
most favorable conditions for this spectral interval. The radiometer spectral response 
(fig. 6(b)) must be considered when discussing the profiles shown in figure 9 because 
atmospheric absorption bands exist at each end of the spectral region. A carbon dioxide 
absorption band centered at 15 pm extends to about 12 pm. Energy from this band, along 
with radiation from ozone in two of its bands at 9.6 pm and 14.08 pm, accounts for the 
slow increase in radiance at tangent heights from 50 to 15 km. This feature of the profile 
is evident in both the measured data Bind the theoretical profile calculated for the condi- 
tions existing during the data period. When the radiometer line of sight sweeps to a tan- 
gent height of about 5 km, the instrument receives radiation from dust and haze present 
at low altitudes. When the field of view contains the earth, blackbody radiation from it is 
also received. 

The relatively smooth increase in all radiance profiles below a tangent height of 
0 km indicates that no clouds came within the field of view of the radiometer. This result 
agrees with the nephanalysis shown in figure 10 which showed few clouds over the region 
of coverage. 

According to the profiles, the 10 to 14 pm interval appears to be well suited for 
horizon sensing because the gradient between 0 km and 10 km is steep and reasonably 
invariant. However, in this interval the presence of clouds could cause extreme varia- 
tions in radiance level. Hence, while it can be concluded that the 10 to 14 pm spectral 
interval offers excellent possibilities for use of horizon sensors in clear weather, it has 
been shown analytically (ref. 2) that the effects of low- altitude weather variations make it 
of doubtful value for general use. 


Flight 3 

For a far infrared sensor whose spectral response is contained entirely within a 
strong atmospheric absorption band, the radiant energy received is primarily thermal 
radiation from the atmospheric gas. In practice, a combination of thermal radiation from 
the atmospheric gas and from the earth or clouds below it is usually obtained because the 
gas is not totally opaque. The selection of a spectral interval which contains a strong 
atmospheric absorption band results in shielding the earth and cloud surfaces by the emis- 
sion from the high-altitude atmospheric gases and, thus, in minimizing the effects of 
variable cloud conditions. 
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Figure 9.- The 10 to 14 ixm horizon radiance profile. Flight 2; May 28, 1%5; observer altitude, 62 km. 

Carbon dioxide, water vapor, and ozone have absorption bands in the infrared 
spectrum and radiate electromagnetic energy in the spectral intervals where they absorb. 
(See refs. 5 and 7 for a more complete discussion of the radiation characteristics of the 
atmosphere.) Of course, the thermal radiation depends on the temperature of the gas in 
question and, consequently, on altitude as indicated in figure 7. For a given pointing direc- 
tion of the radiometer, the optical path from which energy is received can be determined. 
The effective mass of absorbing gas in the path can be found if the pressure and mixing 
ratios of the gases are known. The power radiated by the gas in the spectral interval of 
the sensor can then be determined from Planck's equation if the temperature of the gas 
is also known. The total energy received by the sensor is found by summing the gas 
radiation from each portion of the optical path, while correcting for the attenuation by 
absorption between the radiating point and the sensor. If the optical path is intercepted 
by an opaque body such as a cloud or the earth, the attentuated radiation from the body 
must also be included. A detailed discussion of this procedure is given in reference 2. 

Figure 11 shows the radiance profile obtained in the 14 to 20 pm interval. The 
solid curve is the averaged radiance profile as described for figure 7. The radiance 
elements in this region between 14 qm and 16 (im are CO 2 and in the neighborhood of 
18 qm to beyond equipment cutoff are rotational H 2 O. 


17 



Hephanalysis symbols; 

Cloud! amount: 

0 *= Open; < 20% cloud cover 
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Figure 10.- Nepharta lysis (cloud chart) for 10 to 14 pm spectral interval. Flight 2 ; 18:00 GMT, May 28, 1965. 










Tangent height r km 

Figure 11.- The 14 to 20 pm horizon radiance profife. Flight % June 29, 1%5; observer altitude, 76 km. 


Because the amount of water vapor in the upper atmosphere beyond 25 km is small, 
most of the radiation in this region is due to CC>2* The CO2 region is thought to be a 
most desirable region because the mixing ratio is constant and the absorptivity is great 
enough to screen out the lower atmospheric variability due to weather phenomena. 

The comparison between the results above a tangent height of 25 km shows good 
agreement between theory and measurement and the boundary of all the measurements 
shows the smallest variation of any spectral region measured. The boundary is within 
the ±4 -km uncertainty due to the position and radiance error of the experiment. The 
gradient in this range is steep and thus indications are that the 14 to 16 fim interval is 
usable for horizon sensors. 

The measurements for tangent heights below 20 km show more variability and some- 
what less agreement with theory. The average water-vapor mixing ratio can be estimated 
for a relatively large area but is harder to predict exactly for any specific locality. This 
uncertainty can account for the wider divergence and also the greater scatter in the mea- 
sured data. In addition, some clouds were indicated by the motion pictures but were not 
observed in the measured radiance profiles. These clouds must have been widely 
scattered as indicated by the nephanalysis of figure 12. This disagreement also indi- 
cates variability in water -vapor mixing ratio on a local level. 
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CONCLUSIONS 


Experimental data of the earth's horizon in three spectral intervals have been 
obtained. The resolution (field of view) of the data at the tangent point is approximately 
2 km, and the spatial position of the horizon profile is known to within 4 km. 

The 14 to 20 jLtm interval was shown to be the best spectral interval for horizon 
sensing because it exhibits the least variance at adequate radiant intensity. The 0.8 to 
2.8 jum and 10 to 14 (im spectral intervals do not appear suitable for accurate horizon 
sensing because they are sensitive to low- altitude meteorological changes. 

The experimental data obtained verified the theory used to predict, within experi- 
mental accuracy, the far infrared radiant energy due to carbon dioxide. The use of theory 
for spectral regions where high-altitude water-vapor data are required is questionable 
because no mixing- ratio information is available. The primary scattering theory cannot 
predict accurately the radiant energy in the near infrared region (< 3.0 pm). 

The X-15 airplane provides a convenient test bed for programs requiring high alti- 
tudes and altitude and position data. Also, good weather data are available for the points 
viewed by the radiometer. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 15, 1968, 

715-02-00-01-23. 
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APPENDIX A 


DERIVATION OF THE EXPRESSION RELATING TANGENT HEIGHT 
TO MEASURED PARAMETERS 

In this section the equation expressing tangent height as a function of known flight 
parameters will be developed. It is helpful to define four primary coordinate systems: 
a space -fixed coordinate system, an aircraft coordinate system, a radiometer coordinate 
system, and a local earth coordinate system. The space-fixed coordinate is defined to 
have an x s -axis positive south, a y s -axis positive west, and a z s -axis positive downward. 
The aircraft coordinate system, which is 
assumed to be initially alined with the space- 
fixed coordinate system, has an x a -axis posi- 
tive in the direction of the aircraft nose, a 
y a -axis positive along the right wing, and a 
z a -axis positive downward. The radiometer 
has its x R -axis in the scan plane of the radi- 
ometer line of sight, its y R ~axis in the direc- 
tion of the right wing of the X-15 and lying in 
the plane of the aircraft x a and y a body 
axes, and a z R -axis completing the right-hand 
triad. The local earth coordinate system con- 
sists of an x-axis positive south, a y-axis 
positive east, and a z-axis directed along the 
local vertical, positive upward. 

The scan angle of the radiometer, the 
pitch mounting angle, and the yaw mounting 
angle must be considered to be transformations, as shown in figure 14. In this figure t^ m 
is a positive Euler rotation in the direction shown. The a rotation is a positive Euler 
angle in the direction shown and a m is a negative angle as described. The line-of-sight 
vector, a unit vector initially along the x R -axis, can be described by 



Figure 14.- Radiometer orientation in X-15, 
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The transformation matrices which define the rotations \)/ m , a , and a m are then 
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The alinement of the radiometer coordinates with the aircraft coordinates can be accom- 
plished by the following transformation matrix: 
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Figure 15.- Order of relation through 
attitude angles. 


The angular relationships between aircraft 
coordinates and space coordinates are the stable plat- 
form parameters, roll, pitch, and yaw. The ordered 
rotations indicated in equation (A 5) are used to define 
aircraft coordinates in terms of space coordinates. 
However, it is desired to define space coordinates in 
terms of aircraft coordinates so that the inverse rela- 
tionships must be used. If yaw, pitch, and roll are 
designated by i^, 6, and </>, respectively, in fig- 

ure 15, and if the directions of rotation are positive 
as shown and represent the ordered rotation neces- 
sary to transform space coordinates into aircraft 
coordinates, the space coordinates are found by the 
following transformation: 
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where 


V 


x a 


fx 1 
x a 

Ys 

■ MM>] 

y a 

= m 3 

Ya 
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z a 


z a 


01 - 


cos -sin ip 0 
sin \p cos ip 0 
0 0 1 


0> 


cos 6 0 

0 1 


sin 0 
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|^-sin 6 0 cos ej 


(A6) 


(A7) 


(A8) 


0> 


1 0 

0 COS <P 


0 

-sin <p 


0 sin (j) cos <[> 
Substitution of equation (A4) into equation (A6) yields 


s 

y s 

z a 


= [ M 3][ M 2][ M l] 


(A9) 


(A10) 


The single transformation which will cause the space-fixed and local earth coordinate 
systems to become coincident is represented by 

"l 0 o” 

Mi = | 0 -1 0 (All) 

_0 0 -1_ 

Thus, the total transformation which is obtained by combining equations (A10) and (All) is 

.1 


* MMMM 


(A12) 
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Multiplying matrices in equation (A12) yields the direction of the radiometer line of sight 
in local coordinates as follows: 

cos j//j~cos 8 cos cos^o; + a m ) + sin ejsin <f> sin cos (a + 

+ cos 4> sin{a+ " sin ^[~ cos <p sin cos^y+ 

- sin 4> sin + a m)l 


X R 

Z R 


(A 13) 


-sin ii/j - cos 8 cos \j/ m cos^o;+ + sin ejsin t p sin ipm eos^a + « m J 
+ cos 4> sin^o+ a m)jJ ~ cos ^jcos <p sin s^ m cos^a-f a m ^ 

- sin <p sin(ot + ar m Jj 

- sin 8 cos \f/ m cos (a + - cos s|sin $ sin cos^y + 

+ cos 4> sin(a+ oi m jj 

The direction cosine of the line of sight (fig. 1) is given by 

cos /3 = (-z component of line of sight) 
or 

cos fi = sin 9 cos cos (a + a m ) + cos ej^sin <j> sin cos (a + a m ) + cos <f> sin {a + a m J 

(A14) 

{A 15) 


From figure 1, it follows that 


sin (8 - 


h + r 


h X-15 + r 


Solving for h gives 


h = ^h x _ 15 + rjsin 0 - r 


(A16) 
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where 0 is obtained from equation (A 14). Substituting the expression for 0 into equa- 
tion (A 16) gives 


h = 


sin^cos -1 jsin B cos i^ m cos (a + a m ) + cos s|sin cj> sin ip m cos (a + a m J 


+ cos © sin 


in(a! + a ra )jjj 


( h X- 15 + r ) " r 


(A 17) 
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RADIOMETER DESIGN AND CALIBRATION 


This appendix gives a more detailed description of the radiometer design param- 
eters and of a method used to correct the measured data for the radiometer frequency 
response, as well as a description of the radiometer calibration. 

The smallest obtainable field of view, which determines maximum data resolution, 
is dictated by the resolution limit of the optics and the detectivity of available detectors. 
By using Rayleigh's criterion, for the 12.7-cm aperture of the X-15 radiometer and the 
shortest wavelength considered (X ~ 1.0 j urn), the diffraction- limited resolution is 

% = 1.22 ^ x 57.3 = 5.52 X 10' 4 deg (Bl) 

Because it is almost impossible to produce diffraction-limited optics with an f-number 
of 1, it must be expected that aberrations degrade this resolution. In addition, the flat 
optical filter in the converging beam further degrades the resolution, inasmuch as rays 
incident at different angles from the normal will come to focus at slightly varying focal 
lengths. Resolution tests of the flight optics using visible light at finite and infinite dis- 
tances yielded a telescope resolution estimate of 0.03°. The detector size was selected 
to yield a field of view of 0.13° by 0.13° at the half-response points so that resolution is 
largely determined by the detector. The spectral response of the radiometric systems as 
a function of wavelength for each mission is shown in figure 6. The optical filter and 
detector response were measured separately, multiplied, and corrected for scan and 
objective mirror reflectance to yield the curves in figure 6, 


Detector and Electronics 


The detector subsystems of the radiometer for each mission are somewhat different. 
The characteristics of the three systems are given in table I; however, the parameters 
for only one system (flight 1) will be derived to illustrate the approach used to generate 
the table. The equipment used in flight 1 operated in the 0.8 to 2.8 jam interval and the 
detector was an uncooled lead sulfide photoeonductive cell. Its response combined with 
that of the optical filter gives the overall radiometer response as shown in figure 6, The 
noise equivalent power (NEP) of the detector, in w/hz"' 7 ^, or the root- mean- square value 
of the radiant energy flux, in watts, necessary to give a signal- to- noise ratio of unity, is 


HAV n 

defined as (ref. 7) ^ 


For the 0.8 to 2.8 (jtm interval at a source temperature of 


500° K, at a chopping speed of 15 Hz, and with an equivalent noise bandwidth of 1 Hz, the 
NEP was measured to be 2.07 x 10" ^ w/Hz*/^. The electrical bandpass of the 
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Figure 16.- Transfer function of 0.8 to 2.8 pm radiometer electronics. 


radiometer system electronics is shown in figure 16. The equivalent noise bandwidth is 
defined as (ref. 8) 

B = ~if l G M 2di t B2 i 

1C * 0 


where G(jto) is the transfer function of the system and K is the midband gain of G(jw). 
For the transfer function shown in figure 16, 

B = — 220 = 130 Hz (B3) 

16 

For T * 500° K, f c = 15 Hz, and B = 130 Hz, the NEP is equal to 2.37 x 10" 10 w/hz 1 / 2 . 
The quality of the radiometer may be expressed independently of area by its detectivity, 
given as 


D* = 


^(detector) 


1/2 


Noise equivalent power 


(B4) 


The noise bandwidth is always taken to be 1 Hz. Then the detectivity for the near infrared 
radiometer is 1.45 x 10® cm-HzV^/w. 

The detector is connected in a bridge circuit as shown in figure 17. The active cell 
is exposed to the radiant flux received by the radiometer, while the compensator cell is 
completely shielded from radiant energy. Two cells are selected whose electrical 
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Figure 17.- Detector bridge circuit. 


To amplifier 


characteristics are matched, and both are mounted in close thermal proximity. Gross dc 
unbalance due to changes in ambient temperature is thus eliminated because both legs of 
the bridge change by the same amount. However, because perfectly matched detectors 
are not available, ambient- temperature variations during flight are sufficient to result in 
dc offsets comparable to the radiometric signal output. To filter out these unwanted low- 
frequency signals, the amplifier uses ac coupling, with the low-frequency breakpoint at 
0.01 Hz. A second advantage of using ac coupling is noise filtering, sinee a large part of 
the noise generated in the detector is current noise, which has an intensity variation 
roughly inversely proportional to frequency. The limiting consideration on the choice of 
low-frequency cutoff is dictated by the amount of distortion to the data due to excessive 
droop. 

Because the primary goal of the program is to measure and compare theoretical 
and experimental radiance profiles, the effect of the electronic system must be accounted 
for in the measured radiance. The method for estimating the effect of ac coupling and 
the equivalent electronic- system time delay is discussed in the following paragraphs. 

Figure 18 which gives an assumed trapezoidal radiance input to the radiometer as a 
function of time shows the effect of amplifier time response. Also shown is the radiom- 
eter output waveform for the assumed radiance input. Because of the electronic frequency 
response the output lags behind the input, the output signal does not reach full amplitude, 
and for constant input the output decays to zero. The severity of these factors is deter- 
mined by the low- and high-frequency cutoffs of the system. During data reduction the 
amount of droop that occurs during a horizon crossing must be removed from the record. 
Because the input forcing function is zero while the radiometer looks at space, an expo- 
nential decay of the amplifier output results. Because superposition holds, the exponential 
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Space Horizon Earth 



can be extrapolated through the time of the horizon crossing and its effect removed from 
the recorded signal. Therefore, the correction can be made especially well for those 
profiles where the radiometer line of sight sweeps from space to earth. 

The following procedure can be used to allow for extrapolation of the exponential. 
The amplitude f(f ^ and the time tj - 43 can be measured as indicated in fig- 

ure 18, Because the amplifier has a single low-frequency time constant the output func- 
tions for space input at tj and t 2 are 



(B5) 

-to It 


f (t2) = v i e ' 

(B6) 

t l/ T _t l/ T 

Multiplying equation (B6) by e ' e ' gives 


1 (t 2 ) - v/^/Y* 1 / 7 

(B7) 

It 


f(t 2 ) = Kt(t!) 

(B8) 
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Hence, 

%) “ ffe) = f(*l) - Kf(tx) = f(ti)(l - K) (B9) 


Dividing by (1 - K) gives 



f (‘l) - f (‘2) 
1 - K 


(BIO) 


Equation (BIO) expresses f^t as the difference of the function at two known times, t j 
and t 2 - This result was used to calculate successive points of the exponential* The 
exponential was then subtracted from the flight record to obtain the corrected crossing. 

The time constant of the detector and the location of the upper half-power frequency 
point of the amplifier cause the output to be delayed in time with respect to the input, as 
well as distorted at the beginning and at the end of the ramp (ref* 9), The delay can be 
determined from the detector time constant and the amplifier upper half-power frequency 
point* For an n-pole system, the delay is given as 


At 


n 

2rf c 


(Bll) 


where the corner frequency f c is the half-power frequency point of each pole* Hence, 
for the 3-pole system employed in the X-15 radiometer 


^delay 


= Detector time constant + 


3 

277f e 


(B 12) 


or 


^delay 


- 4*94 msec 


(BIS) 


Once the amount of delay is known it can be subtracted during data reduction so that the 
correct relation is reestablished between the radiometer line of sight and the radiance 
measurements* This correction enables the altitude of the horizon profile to be fixed with 
respect to earth. As indicated previously the characteristics for the far infrared flights 
are similar to those of the near infrared mission* Parameters for the far infrared flights 
are shown in table I, The sensitivity of the thermistor bolometer is two orders of magni- 
tude below that of lead sulfide, and because of the higher low-frequency noise the low- 
frequency cutoff of the amplifier was increased by a factor of 20, 


Calibration 

The exact voltage output of the radiometers as a function of radiant power incident 
on the entrance aperture must be determined* The responsitivity of the radiometer is 
affected by such factors as mirror reflectivity, stray radiation, minute misalinement of 
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the optical elements, and other parameters which are difficult to control and calculate. 
Therefore, the radiometers must be calibrated when assembled as a unit. Three basic 
methods using blackbody sources of known spectral radiance R wX were used to calibrate 
the radiometers. These blackbody reference sources have cavity temperatures control- 
lable and known to ±1° C or better, and have emissivities of 0.98 or greater. Each of the 
three methods was used to calibrate the radiometers, and by cross-checking the data 
obtained from each method confidence is gained in the accuracy of the results. Before 
each method can be discussed, a general equation for radiative power transfer between 
two surfaces must be derived. 

The spectral radiant power emitted from an elemental area dA^ gur j ace i) at tem- 
perature T (g urf a Ce 1) havin & a spectral radiance R wX ( surface into a solid angle 
do? is (ref, 7) 

4 

d P X(surface 1) “ R u>X(surface l) dA (surface l) duJ 


4 

The notation d P x (surface 1) indicates f° ur integrations, two over the area A ( sur f ace 
and two over the solid angle cu. The term R^( sur f ace ^ is the spectral radiance of 

surface 1, which can be a blackbody source. If the temperature and emissivity of the 
blackbody source are known, the spectral radiance is obtainable from either normalized 
blackbody tables or a radiation slide rule. The spectral radiant power emitted by 
dA (surface 1) ^hat * s m1: erce P* ;e cl by a second surface is 


d p x(sur£ace 1 to 2) R cuX(surface l) dA (surface l) dcj (surface 2) 


(B15) 


where do^ gur j ace is the solid angle subtended by the receiving second surface at 
dA (surface 1)' ^ opt ica l filters, such as those in the radiometers, are interposed between 

source and receiver, their spectral transmission jj must be included and equation (B15) 
becomes 


d P X(surface 1 to 2) 77 ^wX(surface l) dA (surface l) duj (surface 2) 


(B16) 


This equation may also be written as 

d p X(rec elver) ~ ^ R wX (source) dA (sour ce) dw (receiver) ^ 

Equation (B17) will now be used to discuss the three calibration methods. 

Collimated-b eam method .- Figure 19(a) shows the experimental setup for the 
eollimated-beam method of calibration. The radiation is parallel (collimated) for the 
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(b! Simplified schematic of radiation transfer. 
Figure 19,“ Collimated-beam method of calibration. 


path between the collimator mirror and the radiometer objective. Therefore, from the 
standpoint of radiation transfer, the configuration may be changed to that of figure 19(b) 
for which equation (B17) is true. The area of the source in equation (B17) is the eff ective 
area of the blackbody source, which includes all points of the source from which radiation 
reaches the detector. By using simple ray tracing the effective area of the blackbody 
source is seen to be the area of the detector multiplied by the magnification of the optics. 
Then by assuming that this image of the detector is smaller than the source, the fol- 
lowing relation is true: 

F c 2 

dA (source) " dA (detector) Z 2 ( B18 ) 

r o 
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Also, ^(receiver) * s so ^ angle subtended by the radiometer objective at the 
source expressed as 

dw, . , = dA (objective) 

(receiver) ^ 2 

*c 


Now, 


d F X(detector) T ^o>X(source) dA (detector) F 2 


F c 2 dA (objective) 


which can be simplified to 


d 4 P 


X(detector) 


_ ^qjX (source) d A (detec tor) dA (obi ec tive) 


(B19) 


(B20) 


(B21) 


By definition, the solid angle da> v of the radiometer is — * [ detl ector ^ therefore 


^ F X(detector) ” ^^■(jX(source) dA (objective) d£i ’ v 
The spectral radiant power on the detector is then expressed as 


(B22) 


F X (detector) _ vR u>\ (source) JJ sur f ace dA (objective) iT so , id d ^ v (B23) 


solid 

angle 


or 


F X (detector) x (source) A (obj ec tive) w v 


(B24) 


The electrical output of the radiometer is measured as a function of radiant power on the 
detector calculated by using equation (B24). This measurement may be repeated for 
several temperatures of blackbody sources so that the radiometer responsivity as a func- 
tion of energy level is determined. 


Extended- area bla ckb o dy- source method .- The extended- area blackbody -source 
method requires a source larger than the objective of the radiometer. The radiometer is 
focused at infinity; therefore, the objective projects each elemental area of the detector 
in a parallel, unchanging beam between objective and source. The dimension l in fig- 
ure 20 can, therefore, be made zero without changing the problem. Then in equation (B17), 
dA( S0U r C e) ls ec l ual to ^(objective)* ^ so ’ the size of the receiver as seen from the 

source dw, ece iver) is (detector) ^ Equations (B21) and (B24) are repeated here for 

F o 

convenience: 
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d P X(detector) 


^ R fc)X(source) d ^(defcector) d ^(objectiv e) 


Po 2 


and 

P X(detector) “ 77 ‘ R, wX{source)'^'(objective) £x,v 

As discussed for the eollimated-beam method, radiometer responsivity as a function of 
irradiance can be found by using equation (B24). This method of calibration is simple to 
implement in the laboratory because the radiant power on the detector is independent of l. 
A problem with this method is that it is very difficult to obtain an extended- area blackbody 
source with constant temperature and emissivity over its entire radiating surface. 



Radiometer 

(simplified) 


Figure 20.- Extended-area blackbody- source method of calibration. 


The Jones method.- The Jones method may be considered to be a special case of the 
extended-area blackbody- source method. As before, optical ray tracing can be used to 
show that in calculating radiant energy transfer, dimension l in figure 21 may be set 

equal to zero. Also dw ( rece i ver ) = — ( det Q ctQr ) > and equation (B17) yields 


An _ 77R c t >X(source) dA (source) dA (detector) 

F X(detector) p 2 


dA 


where — (detector) _ Therefore, 


F o 2 


and 


rk 

d P X(detector) ” ^ R o)X(source) da,v dA (source) 
P X(detector) “ T?R wX(source) aJvA (source) 


(B25) 


(B26) 


(B27) 


36 




ft I 


1 1 mu i 


mi III 


! mm ii 



APPENDIX B 



The advantage of this method is that use can be made of a small -diameter blackbody 
source which is more readily available than extended-area sources. As a result of the 
geometry, only a small portion of the collecting area of the optics is used, and, thus, the 
spectral radiant emittance % of the blackbody must be higher than in the previous 
methods to obtain the same irradiance on the detector. Because of the high-energy den- 
sity of the beam the effects of stray radiation may be appreciable and thus reduce the 
accuracy of the method. The results of the calibration of the radiometers are shown in 
figure 22, The data from the 10 to 14 pm and 14 to 20 pm spectral intervals were sup- 
plied to Langley Research Center by Lockheed Missiles & Space Company under NASA 
Contract NAS1-4932. The power on the detector, in watts, can be obtained from the effec- 
tive radiance used in figures 7, 9, and 11 by multiplying the effective radiance by the 
aperture area and the field of view of the radiometer. The results of the collimated- 
beam method and the extended-area blackbody- source method are seen to agree quite well 
for the three spectral intervals. The results of the Jones method are shown to be in dis- 
agreement, which can be explained by the difficulties with this method previously men- 
tioned, For interpretation of the flight measurements, calibration results from only the 
collimated-beam and extended- area blackbody- source methods were used. 
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Figure 22 .- Radiometer responsitivity for each spectral interval. 
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METEOROLOGICAL DATA AND METHODS USED IN 
DERIVATION OF MODEL ATMOSPHERES 

Meteorological data are required for the evaluation of the horizon data of the radi- 
ometer. These data include temperature and water -vapor mixing- ratio profiles from the 
surface to 50 km or higher and cloud-cover charts (nephanalyses) for the general area 
viewed by the radiometer. A discussion of the data sources, the raw data, the tempera- 
ture and water-vapor vertical profiles, and essential nephanalyses are given in this 
appendix. 


Data Sources 

Te mperature .- Rawinsonde data from U.S. Weather Bureau and military weather 
stations in the X-15 flight area are available daily for heights up to approximately 30 km. 
Meteorological Rocket Network (MRN) soundings are made regularly on Monday, 
Wednesday, and Friday at the White Sands Missile Range, New Mexico, and at Point Mugu, 
California. When the X-15 flight is made on other days, special rocket soundings are 
requested. Thus, temperature and wind measurements for heights up to 50 km or higher 
are made. 

M ixing ratio .- The mixing ratio of CO2 is assumed constant at 314 parts per mil- 
lion, The mixing ratio of water vapor may be derived from the dewpoint data gained 
from raw insondes for heights up to about 10 km. Above this altitude, no humidity data is 
obtained on a regular basis. 

Clouds.- Horizontal weather depiction charts of the U.S, Weather Bureau, and 
nephanalyses from ESSA T s National Environmental Satellite Center served as the bases 
for the modified nephanalyses presented in this appendix. 

Standard surfac e temperature .- The hourly reports from each weather station in 
the scan area were used to derive the standard surface temperature (4 ft above ground 
level). 


General Meteorological Conditions 

The standard meteorological charts showed that for each flight a high-pressure cell 
was dominating the flight area. Therefore, the area was largely free of clouds in each 
case. But knowledge of even scattered cloud conditions is essential, especially in the 
analysis of radiance data for the 0.8 to 2,8 pm and the 10 to 14 pm interval flights. 
Therefore, simplified nephanalyses were prepared with the aid of U.S. Weather Bureau 
horizontal weather depiction charts and Tiros nephanalyses for each flight day. Viewing 
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the motion-picture films from each mission added knowledge of small-scale cloud struc- 
ture to the general picture given by the more standard weather maps. Because the radi- 
ometer in the 10 to 14 pm spectral interval "sees” to the earth's surface, a plot of sur- 
face temperature is included on the nephanalysis for that flight {fig. 10). The results of 
the cloud data for the flights are shown in figures 8, 10, and 12. 

Determination of Temperature and Water-Vapor Profiles 

Figure 23 shows the surface trace of the nominal flight path for the X-15 mission 
and describes the geographical extent of the volume of atmosphere which was sampled on 
a nominal X-15 sensor mission. In constructing a model that is characteristic of the 
atmospheric state during each flight, it is necessary to gather available data describing 
the state from 0 to 50 km and to interpolate these data to describe the volume being 
sampled. Therefore, both rawinsonde and rocketsonde data are needed because the rawin- 
sonde data describes altitudes up to 30 km and the rocketsonde data, those above 30 km. 
Figure 23 shows the rawinsonde sites chosen were Salt Lake City, Utah, and Las Vegas, 
Nevada. The nearest rocketsonde sites are Point Mugu, California, and White Sands 
Missile Range, New Mexico. Because all flights occurred between 12:00 and 00:00 
Greenwich Mean Time (GMT), the derivation of a model atmosphere at rawinsonde alti- 
tudes is amenable to temporal interpolations between 12:00 and 00:00 GMT and spatial 
interpolations between Salt Lake City and White Sands Missile Range. 

The rocketsonde sites are displaced geographically farther from the scan area than 
is desirable. Nevertheless, a spatial interpolation and extrapolation of rocketsonde data 
was performed to yield an estimate of conditions over the scan area at levels higher than 
30 km. Where overlap in rocketsonde and rawinsonde data occurred, the rocketsonde data 
were adjusted to the rawinsonde data. Where the published MRN data lacked a tempera- 
ture correction, the Wagner temperature-correction criterion was applied (refs. 10 
and 11). 

The basic data from both rawinsondes and rocketsondes are plotted in figures 24, 

25, and 26 for flights 1, 2, and 3, respectively. In flight 1, data from only one rocketsonde 
are available. In this flight, of course, no interpolation of rocket data was possible and, 
therefore, the model data derived are less accurate. 

In figures 24, 25, and 26, pressure is plotted against measured temperature. Both 
rawinsonde and rocketsonde data for each flight are plotted on the same chart. It should 
be noted, however, that for the rawinsonde data, pressure is a directly measured quantity, 
whereas for the rocketsonde, it is a derived quantity. The accuracy of all temperature 
data plotted is estimated to be as follows: 

For rawinsonde data 

±1° K up to 12 mbar (nominal 30 km) 
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and for rocketsonde data 

±2° K at 12 mbar (nominal 30 km) 

±3° K at 3 mbar (nominal 40 km) 

±6° K at 0.8 mbar (nominal 50 km) 

±9° K at 0.2 mbar (nominal 60 km) 

For similar studies, see references 10 and 11. 

Calculation of the Water-Vapor Mixing Ratio 

Because the meteorological rockets carried no moisture- sensing instruments, the 
sole water-vapor data available were obtained from the Las Vegas and Salt Lake City 
rawinsonde data. These gave dewpoint information up to approximately 300- mbar alti- 
tude (nominal 10 km). Above this altitude, the nature of the water-vapor content of the 
stratosphere is unknown, and existing observations have given rise to contradictory inter- 
pretations (ref. 12). Therefore the values derived here must be regarded only as esti- 
mates based on one of the existing models. The water-vapor mixing ratios against alti- 
tude were derived by the following procedure: 

(1) Plots of dewpoint against altitude were made for 00:00 and 12:00 GMT for both 
stations, and a mean curve was then selected. 

(2) The saturation vapor pressure corresponding to this dewpoint is read from the 
Smithsonian Meteorological Tables (ref. 13). Saturation vapor pressure over water was 
used for temperatures i -40° C, the saturation vapor pressure over ice was used for 
temperatures <-40° C. 

(3) The mixing ratio w is given by w = 621 p where e is the water-vapor 
pressure at the dewpoint temperature at the given altitude. The units are grams of 
water vapor per kilogram of air. 

(4) Above the peak altitude of moisture data in each case, assumed curves were used. 
These curves were derived in different ways. In the 0.8 to 2.8 pm and 10 to 14 pm flights, 
insufficient moisture data were obtained for any extrapolation. Therefore, in these cases 
the midlatitude mixing- ratio model of the Handbook of Geophysics was used (ref. 14), For 
the 14 to 20 pm spectral interval the mixing ratio profile was obtained by extrapolating 
above the last data point. This profile is parallel to the previously mentioned midlatitude 
model. In every case, above 10 mbar (nominal 31 km) the mixing ratio is assumed con- 
stant to the top of the model atmosphere. 

The resulting mixing ratios are given in the final model atmospheres for flights 1, 

2, and 3 in tables II, in, and IV, respectively. 
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Figure 24.- Pressure-temperature data used in derivation of model atmosphere lor O.S to 2.8 |im spectral interval. 

Flight l; July 8, 1964. 
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Figure 25.- Pressure- temperature data used in derivation of model atmosphere for 10 to 14 pm spectral interval. 

Flight 2 ; May 2S, 1965. 
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Figure 26. Pressure-temperature data used m derivation of model atmosphere for 14 to 20 fim spectral interval 

Flight 3 : June 29, 1965. 


APPENDIX D 


ACCURACY ANALYSIS OF THE EXPERIMENTAL SYSTEM 

The measured radiance profiles are subject to errors introduced by the accuracy 
of the equipment and to errors that arise in reading the flight records. Random errors 
in tangent height are caused by fluctuations in the X-15 pitch-angle and roll-angle data, 
the radiometer scan- mirror-position data, and the accuracy of X-15 altitude radar data. 
Precise knowledge of tangent height is also affected by bias errors such as the mounting 
angle of the radiometer on the X-15 and stable platform mounting and calibration accu- 
racy. The effect of random error due to fluctuations can be estimated by using statistical 
methods; the bias error, estimated from experience. 

To determine the standard deviation of the X-15 attitude data for the present experi- 
ment, a smooth least- squares error curve is fitted to the pitch-angle and roll-angle data. 
Then 5Xjsj, the difference between the fitted curve and the actual curve at each point N 
from this the standard deviation o can be found by using the following equation: 


N 



where N is the number of samples (ref. 15). By using this procedure and the flight 
record, the standard deviations of the pitch and roll angles were found to be 

Oq = 0 . 111 ° 

o, = 0.111° 

The random error of the scan-mirror angle was determined in the same manner, and 
resulted in a a » 0.097°. From experience and from discussion with the X-15 range 
operator, the radar- measurement accuracy is estimated to have a value of 

°h x _ 15 “°* 12 km ' 

The bias errors cannot be measured but must be estimated. The stable-platform 
bias including installation alinement and recording film shrinkage was estimated to be 
° 0 q = 0.08°. The angle between the X-15 center line and the radiometer is determined 
by using a theodolite. The 1 - e accuracy of this measurement is estimated to be 
0.033°. To find how these errors affect the tangent height h use is made of the error- 
propagation equation 
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which applies to linear error analysis when independence is assumed. The term -p is 
obtained from the equation for h and evaluated at a point where the nominal value of 
each parameter is known. The value of h is given by equation (A17) as 



Taking the partial derivative gives 


where 


and 


3h "(Kj s * n ^ + Kg cos 0)(Kf cos 0 ~ Kg sin 0 )(hx_i 5 + r ) 

9 ? = 7 “ 72 

^ 1 - sin 0 + Kg cos By 


= cos cos(a+ o; m ) 


Kg = sin cjf> sin cos^cn + o; ra ) + cos 4> sin(a+ a m ) 


By using the following nominal values: 
B = -4.1° 


- 1°26' 
a= 12 ° 


(D2) 


(D3) 

(D4) 


a m = -2° 17' 

<p = -0.2° 

h X- 15 sl ° km 
r = 6371 km 

The sensitivity coefficient becomes — = 10.98 km/deg. The other sensitivity coeffi- 

0 f j 

cients in equation (Dl) may be found similarly. For the above parameters they are 

P = 0.34 km/deg 

£70 

P = 10.87 km/deg 


4 
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— = o.98 km/km 
9h X-15 

« 10.87 km/deg 

9<*m 

= 10.98 km/deg 
30 C 


By using the errors and sensitivity coefficients, the total estimated standard deviation 
o h is then found to be 1.87 km. If the distribution is assumed to be normal, then 
95.4 percent of the measured points will be within 2 c of the mean. If this is taken as the 
criterion, the accuracy of the data is ±3.74 km. The estimated accuracy will change as 
the flight conditions change. However, for the X-15 the desirable data are gathered 
during the 30- sec range about peak altitude so that changes in the flight conditions are 
small. Thus, ±3.74 km is a good estimate of the experimental accuracy for all the radi- 
ance profiles obtained. 

An estimate of the accuracy of the radiometric calibration may be obtained from 
the data shown in figure 22. The standard deviation a of the individual calibration 
points may be determined in a manner similar to that performed for the aircraft- attitude 
data. For the extended- area and collimated-beam methods and the 14 to 20 pm interval 
radiometer, 2o= 0,0376 pW. Full-scale signal during flight is approximately 0.3 pW, 

n rym 

so that the percent error is x 100 = 12.5. The accuracy of the flight tape recorder 

U i o 

and data- reduction digitizer is thought to be 2 percent, and, thus, the overall radiometric 
accuracy becomes 14.5 percent. 
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TABLE I.- RADIOMETER CHARACTERISTICS 



Flight 1 - July 8, 1964 

Flight 2 - May 28, 1965 

Flight 3 - June 29, 1965 

Spectral region 

0.8 to 2.8 i Lira 

10 to 14 pm 

14 to 20 pm 

Aperture 

12.7-cm diameter 

12.7-cm diameter 

12.7-cm diameter 

Focal length 

12.7 cm 

12.7 cm 

12.7 cm 

Field of view 

0.130 toy o.i3° 

0.13° by 0.13° 

0.13° by 0.13° 

Scan angle 

30° 

30° 

30° 

Scan speed 

20°/ sec 

20°/ sec 

20°/sec 

Filter 

Color filter (3 mm) 

Interference and absorption (1 mm) 

Interference and absorption (1 mm) 

Radiometer relative 
spectral response 

Figure 6(a) 

Figure 6(b) 

Figure 6(c) 

Peak transmission 
of optics 

0.87 

0.62 

0.29 

Detector 

Lead sulfide 
(0.3 mm by 0.3 mm) 

Thermistor bolometer 
(0.3 mm by 0.3 mm) 

Thermistor bolometer 
(0.3 mm by 0.3 mm) 

Noise equivalent 
power 

2.4 x 1Q“ 10 W/W/ 2 

3.4 X 10~ 8 w/hz 1 / 2 

2.5 x 10" 8 W/Hz 1 / 2 

Detectivity 

1.45 x 109 cm-Hz^/w 

1.01 x 10 7 cm-Hz^/w 

1.39 x 10 7 cm-HzVyw 

Electronics bandpass 

0.01 to 220 Hz 

0.2 to 220 Hz 

0.2 to 220 Hz 

Detector time constant 

2.75 msec 

2.8 msec 

2.4 msec 



TABLE II.- MODEL ATMOSPHERE FOR 0.8 TO 2.8 /jm SPECTRAL INTERVAL 

[00:00 GMT, July 9, 1964] 


Pressure, 

mbar 

Temperature, 

o K 

Mixing ratio, 

w, 

g/kg 

850 

300 

5.9 

700 

287 

4.2 

500 

265 

2.6 

400 

255 

.8 

350 

248 

.7 

300 

237 

.05 

250 

227 

.016 

200 

218 

.016 

150 

210 

.013 

a 118 

204 

.010 

100 

208 

.0098 

70 

215 

.014 

55 

219 

.018 

50 

220 

.020 

30 

224 

.045 

20 

230 

.068 

18 

231 

.090 

10 

238 

.150 

7 

241 


5 

257 


3 

261 


2 

272 


1 

278 


.8 

278 


.6 

279 



a Tropopause. 
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TABLE IV.- MODEL ATMOSPHERE FOR 14 TO 20 nm SPECTRAL INTERVAL 

[l8:00 GMT, June 29, 1965] 


Pressure, 

mbar 

Temperature, 

OK 

Mixing ratio, 
g/kg 

850 

299 

4.8 

700 

284 

3.35 

650 

275 

2.5 

500 

267 

.70 

420 

255 

.45 

400 

251 

.41 

300 

235 

.05 

200 

219 

.016 

121 

209 

.010 

a 100 

206 

.0098 

70 

216 

.014 

50 

218 

.020 

30 

224 

.045 

20 

229 

.068 

15 

232 

.100 

10 

234 

.150 

7 

240 


6 

246 


5 

245 


4 

252 


3 

256 


2 

264 


1 

267 


.8 

273 


.6 

272 


.4 

270 


.2 

262 



^ropopause. 
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TABLE V.- TYPICAL MEASURED PROFILES IN0.8TO2.8 jin i SPECTRAL INTERVAL 
[Flight 1; 20:05 GMT July 8, 1964; observer altitude, 50 km] 
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TABLE VII.- TYPICAL MEASURED PROFILES IN 14 TO 20 fxm SPECTRAL INTERVAL 
[^Flight 3; 17:23 GMT, June 29, 1965; observer altitude, 76 kin] 


Tangent 

height, 

km 

Effective radiance, W/m^-sr, for test - 

1 

2 

3 

4 

5 

-15 

4.52 

4.5 

4.61 

4.69 

4.95 

-10 

4.5 

4.48 

4.55 

4.65 

4.91 

-5 

4.45 

4.45 

4.52 

4.63 

4.85 

0 

4.38 

4.35 

4.38 

4.5 

4.72 

5 

4.21 

4.19 

4.25 

4.42 

4.54 

10 

3.72 

3.89 

3.80 

3.75 

4.25 

15 

2.65 

3.05 

2.81 

3.07 

3.75 

20 

2.15 

2.12 

2.25 

2.12 

2.81 

25 

1.88 

2.34 

2.03 

2.25 

2.35 

30 

1.7 

2.05 

1.78 

1.90 

2.06 

35 

1.4 

1.57 

1.55 

1.55 

1.70 

40 

1.0 

1.2 

1.15 

1.05 

1.2 

45 

.55 

.8 

.65 

.5 

.75 

50 

.15 

.25 

.37 

.15 

.2 

55 

0 

0 

0 

0 

0 


NASA-Langley, 1966 21 L-5583 
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